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Abstract 

Discrete symmetries of dynamical flows give rise to relations between pe- 
riodic orbits, reduce the dynamics to a fundamental domain, and lead to fac- 
torizations of zeta functions. These factorizations in turn reduce the labor 
and improve the convergence of cycle expansions for classical and quantum 
spectra associated with the flow. In this paper the general formalism is devel- 
oped, with the A-disk pinball model used as a concrete example and a series 
of physically interesting cases worked out in detail. 
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The periodic orbit theory of classical chaotic dynamical systems has a long 
and distinguished history; initiated by Poincare[[J, and developed as a mathemat- 
ical theory of hyperbolic dynamical systems by Smale, Sinai, Bowen, Ruelle and 
others[|], f|, it has in recent years been applied to many systems of physical 
interest 0, [7], ^ [|. The periodic orbit theory of quantum mechanical systems largely 
parallels this development; originating in the work of Hadamardjl0| and SelbergfTTf 



it has been developed as a quantum mechanical theory by Gutzwiller, Balian and 
Bloch, Berry and others[T2|, |T3 , [TJ|, 15], [15], and has been in the focus of much re- 
cent research. In both the classical and the quantum contexts, one is interested in 
computing spectra of certain evolution operators; this can be done by determining 
zeros of Fredholm determinants or associated zeta functions^, [TTJ, [T7| . The periodic 
orbits emerge in this context essentially through the identity log det = tr log which 
relates the spectrum to the traces of the evolution operators, i.e., the periodic orbits 
or cycles. 

As demonstrated in a series of papers [3,fl8fl-|p2|l, cycle expansions of zeta func- 
tions can be profitably used for the calculation of such spectra in chaotic dynamical 
systems. These systems often come equipped with discrete symmetries, such as 
the reflection and the rotation symmetries of various potentials. We shall show 
here that such symmetries simplify and improve the cycle expansions in a rather 
beautiful way; they can be exploited to relate classes of periodic orbits and fac- 
torize zeta functions, not only in quantum mechanics (where the utility of discrete 
symmetry factorizations is well known |23]), but also in classical mechanics. The 
group-theoretic factorizations of zeta functions that we develop here were first in- 
troduced and applied in ref. M. They are closely related to the symmetrizations 



introduced by Gutzwiller [24| in the context of the semiclassical periodic orbit trace 



formulas, recently put into more general group-theoretic context by Robbins||25||. 



whose exposition, together with Lauritzen'spq] treatment of the boundary orbits, 
has influenced the presentation given here. A related group-theoretic decomposition 
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in context of hyperbolic billiards was utilized in ref. [27 



Invariance of a system under symmetries means that the symmetry image of a 
cycle is again a cycle, with the same weight. The new orbit may be topologically 
distinct (in which case it contributes to the multiplicity of the cycle) or it may be 
the same cycle, shifted in time. In the latter case, the cycle can be subdivided into 
segments, each of which is a symmetry image of an irreducible segment. The period 
or the action of the full orbit is given by the sum along the segments, whereas the 
stability is given by the product of the stability matrices of the individual segments. 
The phase space can be completely tiled by a fundamental domain and its symmetry 
images ( "phase space" in this paper stands for the coordinates of any <i-dimensional 
dynamical system whose evolution is described by a set of first order differential 
equations or iterative mappings, not necessarily Hamiltonian) . The irreducible 
segments of cycles in the full space, folded back into the fundamental domain, are 
closed orbits in the reduced space. 

The main point of this paper is that if the dynamics possesses a discrete sym- 
metry, the contribution of a cycle p of multiplicity m p to a dynamical zeta function 
factorizes into a product over the <i a -dimensional irreducible representations D a of 
the symmetry group, 

(i - t P r p = n det (! - D a (h) p t p ) d - , t p = ti /mp , (i) 

a 

where tp is the cycle weight evaluated on the fundamental domain, g is the di- 
mension of the group, hp is the group element relating the fundamental domain 
cycle p to a segment of the full space cycle p, and m p is the multiplicity of the p 
cycle. Emergence of symmetrized subspaces, a common phenomenon in quantum 
mechanics, is perhaps surprising in a classical dynamics context. The basic idea 
is simple: in classical dynamics, just as in quantum mechanics, the symmetrized 
subspaces can be probed by linear operators (observables) of different symmetries. 
If a linear operator commutes with the symmetry, it can be block-diagonalized and 
the associated determinants will therefore factorize. 

This paper is meant to serve as a detailed guide to computation of zeta functions 
for systems with discrete symmetries. We develop here the cycle expansions needed 
for evaluation of factorized zeta functions, and examplify them by working out a 
series of cases of physical interest: C*2, C^, C% v and C± v symmetries. For instance, 



one has a C 2 symmetry in the Lorenz system Eq, 29], the Ising model, and in the 
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3- dimensional anisotropic Kepler potential [E4], ISO, |3T|], a symmetry in Henon- 



Heiles type potentials |32], [33], ^J, |2q| , a symmetry in quartic oscillators^, [36 
in the pure x 2 y 2 potential [37. 38] and in hydrogen in a magnetic field |]59|, and a 
C*2i; = C 2 x C 2 symmetry in the stadium billiard|25|]. 

We will illustrate our results using the pinball scattering by three and four 
disks |4(| as an example. Besides their intrinsic interest as examples of classical 



and quantum mechanical chaotic dynamics p3|, [41], [4^], they are also relevant to 



smooth potentials. The pinball model may be thought of as the infinite potential 
wall limit of a smooth potential, and, much as the 1-d tent map captures the 
topology of a general unimodal map, the iV-disk symbol dynamics can serve as a 
covering symbolic dynamics in smooth potentials. One may either define potential 
wall collisions in phase space [^, [E|] or one may start with the infinite wall limit and 
continuously relax an unstable cycle onto the corresponding one for the potential 
under investigation. If things go well, the cycle will remain unstable and isolated, 
no new orbits (unaccounted for by the iV-disk symbolic dynamics) will be born, and 
the lost orbits will be accounted for by a set of pruning rules. For example, this 
adiabatic approach has been profitably used in ref. jH| in disproving the conjecture 
that the x 2 y 2 potential is ergodic. Its validity has to be checked carefully in each 
application, as things can easily go wrong; for example, near bifurcations the same 
naive symbol string assignments can refer to a whole island of distinct periodic 
orbits. 

In addition to the symmetries exploited here, time reversal symmetry and a 
variety of other non-trivial discrete symmetries can induce further relations among 
orbits; we shall point out several of examples of cycle degeneracies under time 
reversal. We do not know whether such symmetries can be exploited for further 
improvements of cycle expansions. 

The paper is organized as follows. In the next section, we recall some basic facts 
of the zeta-function formalism and cycle expansions, and describe the symbolic 
dynamics of the iV-disk model. In section |^ we illustrate the utility of zeta function 
formalism by using it to count cycles. In section [| we introduce discrete symmetries 
and apply them to identify degenerate classes of cycles. In section ^] we describe 
the reduction of the iV-disk dynamics to the fundamental domain and the special 
treatment required by boundary orbits. Finally, in section ^] we apply symmetries 
to reduce the symbolic dynamics and to factorize zeta functions. Several examples 
are worked out in detail. Cycle expansions for the (symmetry unreduced) 3- and 
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4-disk dynamics are listed in the appendix. We conclude with a summary and an 
outlook for further work. 



2 Preliminaries 

Here and in sect. [3] we review the cycle expansion formalism; the subject proper of 
the paper, group-theoretic factorizations, commences only in sect. f|. The reader 
might profitably start with sect. |], and refer back to the preliminaries of sects. [2] 
and |3| as need arises. 



2.1 Zeta functions 

Transfer operators and the associated zeta functions are treated extensively in the 
literature 0. Here we merely state the results needed for the purposes of this paper, 
in the notation of ref. ]18| . 



The general setting is as follows: given a dynamical system, presented either 
as a d- dimensional map f(x), or as a d+1 dimensional flow (in the latter case, the 
flow can be reduced to a c/-dimensional mapping by means of appropriate Poincare 
sections) , one is interested in time evolution of certain distributions, such as classical 
probability distributions and quantum mechanical wave functions. The effect of the 
dynamics on such distributions is given by linear evolution operators, such as the 
integral kernel used in evaluation of the escape rate from a repeller[^5| described 
by a map x n+1 = f(x n ), 

C{y,x) = 5{y-f{x)) . (2) 

The time dependence of distributions in question is determined by the eigenspec- 
trum and eigenfunctions of evolution operators, and the problem that concerns us 
here is the problem of effective evaluation of such eigenspectra for a given dynamical 
system. 

The eigenvalues A*, of £ are the inverses of the zeros of 

oo 

Z(z) = Det(l - zC) = H(l - z\ k ) . (3) 

k=0 
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The periodic orbit approach to the determination of zeros of Z(z) is based on the 
observation that the determinant of an operator is related to its traces by 



oo z n 



det(l-z£)=exp tr£ n (4) 



n=l n 



Traces of an operator like @) receive contributions from all prime cycles p of period 
n p and their multiple traversals, weighted by the cycle Jacobians 



^)=exp -EE: ,, T. ,sA ■ (5) 




This expression for Det(l — zC) is specific to the mapping (Q). In more general 
settings [0], z Up is replaced by a weight whose precise form depends on the particular 
average being computed; for example, in the corresponding determinant for classical 
smooth flows, z np is replaced by e sTp , where T p is the p cycle period ^6] . The 



quantum mechanical kernel corresponding to (|2|) is smeared out by a path integral, 
but in the semi classical approximation [14] the zeta function] VJJ has essentially the 



same form as (pj): 

1 e iS p {E)r-i-K^ p r/2 



P r =i r | det (l - 3$) I 1 2 



S p denotes the classical action and \x v is the Maslov index. As the group-theoretic 
factorizations that we shall develop here rely only on the linearity of evolution 
operators, they will apply to both the classical and the quantum cases. 

For evaluation of spectra, the expansion (EH can be used pretty much as it 



standspOj) ft can be expanded as a multinomial in cycle weights^], or the leading 



eigenvalues can be extracted from the associated dynamical zeta function || 

1/C = II( 1 -U tp = z np /A p or t p = e is P (E)-i^/2 /A i/2 ? (g) 

V 

obtained by approximating det (1 — 3 P ) = nf=i (1 — Ap><») by the product of ex- 
panding eigenvalues A p = [fla^al i n (§)■ For example, for a 1-dimensional ex- 
panding map each periodic point contributes with the weight 1/|1 — AT|: 

co n oo -i n p r 

E-^ n = EE- , — * — r • 

n=l p r=l 'I 1 iv pl 
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Substituting 1/|1 - A| = |A| _1 (1 - A" 1 )" 1 = (A^EA^ we obtain 




The dynamical zeta function (||) is the k — part of the above Selberg-type prod- 
uct. The full zeta functions (|5|) are infinite products over dynamical zeta functions 
which depend on both the expanding and the contracting eigenvalues. Most of the 
developments below are independent of the precise form of the zeta functions. 

As the dynamical zeta functions have particularly simple cycle expansions, a 
simple geometrical shadowing interpretation of their convergence, and as they suffice 
for determination of leading eigenvalues, we shall concentrate in this paper on their 
factorizations; the full Z(z) determinants can be factorized by the same techniques. 
To emphasize the group theoretic structure of zeta functions, we shall combine all 
the non-group-theory dependence of a p-cycle into a cycle weight t p . We shall also 
often absorb z into the transfer operator: zC — > C, z Up t p — > t p . 

The first prerequisite for converting expressions like into cycle expansions is 
efficient enumeration of periodic orbits - the problem to which we turn next. 

2.2 Symbolic dynamics 

The key to a theory of a chaotic dynamical system is its qualitative, topological 
description 0, or, as it is usually called, its symbolic dynamics. The strategy is to 
partition phase space into topologically distinct regions, associate with each region 
a symbol from an alphabet, and use those symbols to label every possible trajectory. 
Covering symbolic dynamics assigns a distinct label to each distinct trajectory, 
though there might be symbol sequences which are not realized as trajectories. If 
all possible symbol sequences can be realized as physical trajectories, the symbolic 
dynamics is called complete; if some sequences are not allowed, the symbolic dy- 
namics is pruned (the word is suggested by "pruning" of branches corresponding 
to forbidden sequences for symbolic dynamics organized hierarchically into a tree 
structure). In that case the alphabet must be supplemented by a grammar, a set of 
pruning rules. A periodic symbol string corresponds to a periodic orbit or a cycle. 
Periodic orbits will here be distinguished by a bar over the primitive symbol block, 
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but we often omit the bar if it is clear from the context that we are dealing with a 
periodic orbit. 

These concepts are easily illustrated by the pinball models that we shall study 
here. Consider the motion of a point particle in a plane with N elastically reflecting 
convex disks. Any trajectory can be labelled by the sequence of disk bounces, and 
a symbolic dynamics is given by the alphabet of N symbols {1, 2, 3, ■ • • , N}. As the 
bodies are convex, there can be no two consecutive reflections off the same disk, 
hence the first rule of the grammar of allowed sequences is that symbol repetitions 
_11_, _22_, • • •, _NN_ are forbidden. More generally, we shall refer to a symbolic 
dynamics as an "iV-disk" symbolic dynamics if the phase space can be partitioned 
in N distinct regions such that an orbit starting in a partition can in one step reach 
all other partitions except itself. 

A finite length scattering trajectory is not uniquely specified by its (finite) sym- 
bol sequence, but an unstable cycle (consisting of infinitely many repetitions of a 
prime building block) is. We shall show in sect. ^| that the prime cycles for such 
simple grammars are easily enumerated; for example, table 1 contains a list of all 
prime 3-disk cycles up to length 9, and table 2 contains a list of prime 4-disk cycles. 

An important effect of a discrete symmetry is that it tesselates the phase space 
into copies of a fundamental domain, and thus induces a natural partition of phase 
space. The group elements g = {a, b, ■ ■ ■ ,d} which map the fundamental domain 
M into its copies gM, can double in function as letters of a symbolic dynamics al- 
phabet. If the dynamics is symmetric under interchanges of disks, the absolute disk 
labels = 1, 2, ■ ■ • , N can be replaced by the symmetry- invariant relative disk-^disk 
increments gi, where gi is the discrete group element that maps disk i — \ into disk 
i. Experience shows that more often than not specifics of the model at hand dictate 
the choice of symmetry reduced symbolic dynamics, so rather than attempting to 
develop a general procedure here, we shall demonstrate the reduction for a series 
of specific examples in sect. |6|. An immediate gain arising from symmetry invari- 
ant relabelling is that iV-disk symbolic dynamics becomes (N — l)-nary, with no 
restrictions on the allowed sequences [40, f48|. However, the main gain is in the close 
connection between the symbol string symmetries and the phase space symmetries 
which will aid us in the zeta function factorizations. Once the connection between 
the full space and the reduced space is established, working in the fundamental 
domain (ie., with irreducible segments) is so much simpler that we never use the 
full space orbits in actual computations. 
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Whether this symbolic dynamics is complete (as is the case for sufficiently sepa- 
rated disks), pruned (for example, for touching disks), or only a first coarse graining 
of the topology (as, for example, in systems with islands of stability) depends on 
the details of the dynamics and requires further case-by-case investigation. In the 
iV-disk model outlined above we have tacitly assumed that the disks are sufficiently 
separated so that all possible symbol sequences (excluding symbol repeats) are 
realized as physical trajectories. If the disks shadow each other, further infinite 
families of sequences are pruned [HO, ISO, pTJ. Determining the pruning rules is in 



general a highly non-trivial undertaking, carried out so far for only a few dynamical 
systems 0, [52| . In this paper we assume that the disks are sufficiently separated 
so that there is no additional pruning beyond self bounces. This assumption does 
not affect our main result, the symmetry induced factorizations - it only affects the 
cycle counting of sect. || 

2.3 Cycle expansions 

In the simplest of the cases that we shall discuss here (an example is the fundamental 
domain symmetric subspace of a 3-disk repeller, discussed in sect. |Z. 2|) the system is 
described by a complete binary symbolic dynamics. The Euler product (^) is given 
by@ 

1/C = (I - zt )(l - zh^l - z 2 t 01 )(l - zh 001 )(l - zh ou ) 

(1 — ^ 4 ^oooi)(l ~ z toon){l — ^ 4 ioiii)(i — ^ 5 tooooi)(l ~~ z ^ooon) 

(1 — 2; 5 ^ooioi)(l — £ 5 £oom)(l — £ 5 £oion)(l — £ 5 £oim) • • • (8) 

The cycle expansion is obtained by multiplying out the Euler product and group- 
ing together the terms of the same total string length (same power of z): 

1/C = l-*t -zti-* 2 [(toi-Mo)] 

— z ^ [(^ooi — ^oi^o) — (^on — tmti)} 

— z [(^0001 — ^O^OOl) + (^0111 — ^Oll^l) 

+ (^0011 — ^ooi^i — ^o^on + ^oA)] — ■ ■ • (9) 



The terms grouped in brackets are the curvature corrections |5B|]; the terms 
grouped in parenthesis are combinations of longer orbits and their shorter "shad- 
owing" approximants. In the counting limit t p — 1, and all such shadowing com- 
binations vanish (as we show in more detail in sect. El). The practical utility of 
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cycle expansions, in contrast to direct averages over periodic orbits such as in the 
trace for inula [ 12 1, lies precisely in this organization into nearly cancelling combi- 
nations: cycle expansions are dominated by short cycles, with long cycles giving 
exponentially decaying corrections. 

Further examples of cycle expansions are given in sect. ^ and in the Appendix. 



3 Counting cycles 



In this section we apply the cycle-counting methods of ref. |18|] to the iV-disk prob 



lem. This section is not essential to the main, group-theoretic thrust of this paper, 
but in practice cycle counting is of some use as a check of correctness of various 
cycle expansions. 

There are N n = N n possible distinct strings of length n composed of N letters. 
These N n strings include all M d prime cycles whose period d equals or divides n, 
and each d— cycle contributes its cyclic permutations, d in number: 

N n = J2dM d . (10) 

d\n 



The number of prime cycles follows by Mobius inversion 54 



M n = n- 1 ^^^) N d . (11) 



d\n 



where = 1, fi(n) = if n has a squared factor, and [i{pxP2 ■ ■ -Pk) = (— l) fe if all 
prime factors are different. 

For example, from two symbols {0, 1} one can form M n = 2, 1, 2, 3, 6, 9, 18, . . . 
prime cycles, i.e., there are two fixed points and 1, one prime 2-cycle 01, two 



3-cycles 001 and 011, three 4-cycles 0001, 0011, 0111 (note that e.g. 1010 = 01 2 is 
not prime), and so forth. In table 3 we list the number of prime orbits up to length 
10 for 2-, 3- and 4- letter complete symbolic dynamics. 

For a generic dynamical system not all M n prime periodic symbol strings are 
realized as physical orbits: in such cases M n is only an upper bound to the actual 
number of prime n-cycles. To count correctly, we need to prune the disallowed 
orbits, ie. specify the grammar of the allowed sequences. A simple example of 
pruning is the exclusion of "self-bounces" in the iV-disk pinball. To determine the 
number of periodic orbits, consider a [N x N] transfer matrix whose elements are 
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Ty = 1 if a transition from disk j to disk i is possible, and otherwise. The 
number of points that are mapped back onto themselves after n iterations is given 
by N n = tr T n . For a complete iV-ary dynamics all entries would equal unity 
/l 1 ... 1\ 

; ' - ; • < 12 > 

VI 1 ... 1/ 

and tr T™ = N n . The pruning of self-bounces eliminates the diagonal entries, 
TN-disk —T c — 1, so the number of the iV-disk periodic points is 

N„ 



trn_ disk = (N-ir+(-ir(N-i). 



(13) 



Mobius inversion (|TT|) now yields 



M. 



N-disk 
n 



n 



d\n 
M (N-1) 



AT 



n 



an 



for n > 2 



(14) 



0. The number of periodic points of period 



N, hence there are M?~ disk = N(N - 1) /2 prime cycles of length 2; 



There are no fixed points, Mf~ dtsk 

2 is iV 2 — iV, 11C11UC U11C1C CllC jh 2 

for lengths n > 2, the number of prime cycles is the same as for the complete 
(JV — l)-ary dynamics. 

The simplest application of the cycle expansion of dynamical zeta function @ 
is the evaluation of the topological entropy. The topological entropy is the growth 
rate of the number of orbits as a function of the length of their symbol sequences: 



h = lim In N n /n . 



(15) 



The topological entropy is given by the logarithm of the largest eigenvalue of the 
transition matrix T such as the one given in the above example. Dynamical zeta 
functions enter via the relationship (see eq. (18) of ref. 



det(l-zT) = Y[(l-t p ) 



(16) 



where for the topological entropy the weight assigned to a prime cycle p of length 
n p is t p = z np if the cycle is allowed, or t p = if it is pruned. Expanded in powers 
of z one finds 



c^r 1 = n (! - ^ = i+e^- 

V k=l 



;i7) 
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This function is called the topological zeta functiong, |55fl ; if the grammar is finite, 
it reduces to the topological polynomial. The topological entropy h is given by the 
smallest zero z = e~ h . 

For complete symbolic dynamics of N symbols, the topological polynomial is 
simply 

((zy^l-Nz, (18) 

whence the topological entropy h = In N. 

One consequence of the finitness of topological polynomials is that the contri- 
butions to curvatures at every order are even in number, half with positive and half 
with negative sign. For instance, for complete binary labelling @, 

Q = —^oooi — £0011 — £0111 — totoiti 

+ ^0^001 + ^o^oii + £001^1 + £01 A • (19) 

We see that 2 3 terms contribute to C4, and exactly half of them appear with a 
negative sign - hence if all binary strings are allowed, this term vanishes in the 
counting expression. The number of terms can be counted using the identity 

na+^) = n^- (20) 

p p 1 h v 

Substituting fll8| ) we obtain 

1 iu 7 2 00 

11(1 + **) = T^wr = 1 + Nz + £, fe K - n"- 1 ) . (21) 

p 1 i\z k=2 

The z n coefficient in the above expansion is the number of terms contributing to 
c n , so we find that for complete symbolic dynamics of iV symbols and n > 1, the 
number of terms contributing to c n is (N — l)iV n_1 . 

This technique is easily generalized to cycle expansions for the N-disk symbol 
sequences. Consider for example the 3-disk pinball. The prohibition of repeating a 
symbol affects counting only for the fixed points and the 2-cycles. Everything else 
is the same as counting for a complete binary dynamics (eq fl!H|)). To obtain the 
topological zeta function, just divide out the binary 1- and 2-cycles (1 — zt )(l — 
zt — z 2 toi) and multiply with the correct 3-disk 2-cycles (1 — z 2 ti2)(l — z 2 ti 3 )(l — 

z 2 t 23 y. 

(l-z 2 ) 3 

= (1 -2z)(l + z) 2 = 1 -3z 2 -2z 3 . (22) 
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The 4-disk pinball topological polynomial can be derived in the same way: the 
pruning affects again only the fixed points and the 2-cycles 

(l-z 2 ) 6 

= (l-3z)(l + z) 3 = I - 6z 2 - 8z 3 - 3z 4 , (23) 
and, more generally, for an iV-disk pinball, the topological polynominal is given by 

(1 _ z 2\N{N-l)/2 
1/ (N-disk = (1 - (N - l)z) — z Y~l(l _ z 2yN-l)(N-2)/2 

= (l-(N-l)z)(l + z) If - 1 . (24) 

The topological polynomial has a root z~ x = N — 1, as we already know it should 
from (P^). We shall see in sect. |6] that the other roots reflect the symmetry factor- 
izations of zeta functions. Since the zeta functions reduce to polynomials, we are 
assured that there are just a few fundamental cycles and that all long cycles can 
be grouped into curvature combinations. For example, the fundamental cycles in 
(p2|) are the three 2-cycles which bounce back and forth between two disks and the 
two 3-cycles which visit every disk. It is only after these fundamental cycles have 
been included that a cycle expansion is expected to start converging smoothly, i.e., 
only for n larger than the lengths of the fundamental cycles are the curvatures c n , a 
measure of the deviations between long orbits and their short cycle approximants, 
expected to fall off rapidly with n. 

Conversely, if the grammar is not finite and there is no finite topological poly- 
nomial, there will be no "curvature" expansions, and the convergence will be poor. 



That is the generic case, and one strategy|18, 51, 52] for dealing with it is to find a 



good sequence of approximate but finite grammars; for each approximate grammar 
cycle expansions yield exponentialy accurate eigenvalues, with succesive approxi- 
mate grammars converging toward the desired infinite grammar system. 

The iV-disk topological zeta function may also be used to count the number of 
terms in the curvatures. For example, for the 3-disk pinball we get 

T-r , 1 - 3z A - 2z 6 9 » z 4 (6 + 12z + 2z 2 ) . , 

The coefficients are 1, 0, 3, 2, 6, 12, 20, 48, 84, 184, . . .. That means that, e.g., c 6 
has a total of 20 terms, in agreement with the explicit 3-disk cycle expansion (|62|) 
of the appendix. 
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The above concludes our review of cycle expansions for general dynamical sys- 
tems; now we turn to the main subject of this paper, the role of discrete symmetries 
in cycle expansions. 



4 Discrete symmetries 



A dynamical system is invariant under a symmetry group G = {e, g2, ■ ■ ■ , g g } if 
the equations of motion are invariant under all symmetries g G G. For a map 
Xn+i = f{ x n) an d the evolution operator C(y,x) defined by @ this means 

f(x) = g -1 /(gz) 

C(y,x) = £{gy,gx) . (26) 

Bold face letters for group elements indicate a suitable representation on phase 
space. For example, if a 2-dimensional map has the symmetry X\ — > — X\, X2 — ► — 3?2j 
the symmetry group G consists of the identity and C, a rotation by ir around the 
origin. The map / must then commute with rotations by 7r, f(Cx) = Cf(x), with 
C given by the [2 x 2] matrix 

" -1 

C = . (27) 

-1 

C satisfies C 2 = e and can be used to decompose the phase space into mutually or- 
thogonal symmetric and antisymmetric subspaces by means of projection operators 



Pa* 

£Ai{y,x) 

£>A 2 (y,x) 



;(e + C) 



Pa 



;(e-C) 



P Al £(y, x) = - (C{y, x) + C(-y, x)) , 
P A2 C{y,x) = - (C(y,x) - C(-y,x)) . 



(2f 



More generally[25, 26| the projection operator onto the a irreducible subspace of 
dimension d a is given by P a = (d a /\G\) X«(^)h _1 , where Xa(h) = tr D a (h) are 
the group characters, and the transfer operator £ splits into a sum of inequivalent 
irreducible subspace contributions Satr£ a , 



£ Q {y,x) = j£r Xa(h)C(h 1 y,x) 



(29) 



The prefactor d a in the above reflects the fact that a d a - dimensional representation 



occurs d a times. 
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4.1 Cycle degeneracies 

If g G G is a symmetry of the dynamical problem, the weight of a cycle p and 
the weight of its image under a symmetry transformation g are equal, t gp = t p . 
The number of degenerate cycles (topologically distinct, but mapped into each 
other by symmetry transformations) depends on the cycle symmetries. Associated 
with a given cycle p is a maximal subgroup H p C G, H p = {e,b 2 ,b 3 , . . . ,bh} of 
order h p , whose elements leave p invariant. The elements of the quotient space 
b G G /Ti p generate the degenerate cycles bp, so the multiplicity of a degenerate 
cycle is m p = g/h p . 

Taking into account these degeneracies, the Euler product (|6]) takes the form 

IK 1 -**) = IK 1 -<*)""■ ( 30 ) 

v p 

Here p is one of the m p degenerate cycles, picked to serve as the label for the entire 
class. Our labelling convention is usually lexical, i.e., we label a cycle p by the cycle 
point whose label has the lowest value, and we label a class of degenerate cycles by 
the one with the lowest label p. In what follows we shall drop the hat in p when 
it is clear from the context that we are dealing with symmetry distinct classes of 
cycles. 

4.2 Example: C% v invar iance 

An illustration of the above is afforded by C 3 „, the group of symmetries of a pinball 
with three equal size, equally spaced disks, fig. |l|. The group consists[|6| of the 
identity element e, three reflections across axes {cr 12 , a 2 3, c 13 }, and two rotations 
by 27r/3 and 47r/3 denoted {C3, Cf}, so its dimension is g = 6. On the disk labels 
{1,2,3} these symmetries act as permutations which map cycles into cycles. For 
example, the flip across the symmetry axis going through disk 1 interchanges the 
symbols 2 and 3; it maps the cycle 12123 into 13132, fig. ^|a. 

The subgroups of C^ v are C v , consisting of the identity and any one of the 
reflections, of dimension h = 2, and C3 = {e,C3,Cf}, of dimension ft, = 3, so 
possible cycle multiplicities are g/h = 2, 3 or 6. 

The C3 subgroup invariance is exemplified by the cycles 123 and 132 which are 
invariant under rotations by 27r/3 and 47r/3, but are mapped into each other by 
any reflection, fig. ^b; H p = {e, C 3 , Cf }, and the degeneracy is g/h c& = 2. 
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The C v type of a subgroup is exemplified by the invariances of p = 1213. This 
cycle is invariant under reflection cr 2 3{1213} = 1312 = 1213, so the invariant sub- 
group is Hp = {e, (J23}. Its order is hc v = 2, so the degeneracy is mp = g/hc v = 3; 
the cycles in this class, 1213, 1232 and 1323, are related by 27r/3 rotations, fig. ||]c. 

A cycle of no symmetry, such as 12123, has 7i p = {e} and contributes in all six 
terms (the remaining cycles in the class are 12132, 12313, 12323, 13132 and 13232), 
fig. |a. 

Besides the above discrete symmetries, for Hamiltonian systems cycles may be 
related by time reversal symmetry. An example ||57|| are the cycles 121212313 and 
121212323 = 313212121 which are related by no space symmetry (fig. §1). 

The Euler product (§) for the C$ v symmetric 3-disk problem is given in the 
Appendix, eq. (|62|). 

5 Dynamics in the fundamental domain 

So far we have used the discrete symmetry to effect a reduction in the number of 
independent cycles in cycle expansions. The next step achieves much more: the 
symmetries can be used to restrict all computations to a fundamental domain. We 
show here that to each global cycle p corresponds a fundamental domain cycle p. 
Conversely, each fundamental domain cycle p traces out a segment of the global 
cycle p, with the end point of the cycle p mapped into the irreducible segment of p 
with the group element hp. 

If the dynamics is invariant under a discrete symmetry, the phase space M can 
be completely tiled by the fundamental domain M and its images aM, bM, . . . 
under the action of the symmetry group G = {e, a, b, . . .}, 

M=J2M a = ^aM. 

aeG aGG 

In the above example (|7|) with symmetry group G = {e,C}, the phase space 
M = {x\-X2 plane} can be tiled by a fundamental domain M = {half-plane x\ > 0}, 
and CM = {half-plane X\ < 0}, its image under rotation by tt. 

If the space M is decomposed into g tiles, a function <p{x) over M splits into a 
g-dimensional vector (fi a (x) defined by 4> a (x) = 0(x) if x G M a , (f> a (x) = otherwise. 
Let h = a6 _1 be the symmetry operation that maps the endpoint domain into 
the starting point domain M a , and let D(h)b a , the left regular representation, be 
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the [g x g] matrix whose b, a-th entry equals unity if a = hb and zero otherwise; 
D(h) ba = 5hb, a - Since the symmetries act on phase space as well, the operation h 
enters in two guises: as a [g x g] matrix D[h) which simply permutes the domain 
labels, and as a [d x d] matrix representation h of a discrete symmetry operation 



on the d phase-space coordinates. For instance, in the above example (|27|) h e C 2 
and D(h) can be either the identity or the interchange of the two domain labels, 



Die) 



1 
1 



, D(C) 



1 

1 



(31) 



Note that D(h) is a permutation matrix, mapping a tile M a into a different tile 
Mha 7^ M a if h e. Consequently only -D(e) has diagonal elements, and tr D(h) = 
gSh, e - However, the phase-space transformation h ^ e leaves invariant sets of 
boundary points; for example, under reflection a across a symmetry axis, the axis 
itself remains invariant. The boundary periodic orbits that belong to such point- 
wise invariant sets will require special care in tr£ evaluations. 

One can associate to the evolution operator (Q) a [gxg] matrix evolution operator 
defined by 

C ba (y,x) = D(h) ba C(y,x) , 

if x G M a and y G M b , and zero otherwise. Now we can use the invariance condition 
( |26l) to move the starting point x into the fundamental domain x = ax, £(y, x) = 
£(a _1 y, x), and then use the relation a~ 1 b = h^ 1 to also relate the endpoint y to its 
image in the fundamental domain, C(y, x) = £(h _1 y, x). With this operator which 
is restricted to the fundamental domain, the global dynamics reduces to 

£ba(y,x) = D(h) ba C(y,x) . 

While the global trajectory runs over the full space M, the restricted trajectory is 
brought back into the fundamental domain M any time it crosses into adjoining 
tiles; the two trajectories are related by the symmetry operation h which maps the 
global endpoint into its fundamental domain image. 

Now the traces (§) required for the evaluation of the eigenvalues of the transfer 
operator can be evaluated on the fundamental domain alone 



tr£= / dxC(x,x) = / dx V" tr _D(/i) £(h x x, 



x 



(32) 
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The fundamental domain integral / dx £(h _1 x, x) picks up a contribution from 
every global cycle (for which h = e), but it also picks up contributions from shorter 
segments of global cycles. The permutation matrix D(h) guarantees by the identity 
trD(h) — 0, h ^ e, that only those repeats of the fundamental domain cycles p 
that correspond to complete global cycles p contribute. Compare, for example, the 
contributions of the 12 and cycles of fig. [| tr D(h)C does not get a contribution 
from the cycle, as the symmetry operation that maps the first half of the 12 
into the fundamental domain is a reflection, and tr D(a) = 0. In contrast, a 2 = 
e, tr D(a 2 ) = 6 insures that the repeat of the fundamental domain fixed point 
tr (D(h)C) 2 = 6tg, gives the correct contribution to the global trace tr C 2 = 3 ■ 2ii 2 . 

Let p be the full orbit, p the orbit in the fundamental domain and hp an ele- 
ment of Tip, the symmetry group of p. Restricting the volume integrations to the 
infinitesimal neighborhoods of the cycles p and p, respectively, and performing the 
standard resummationsfl8|l, we obtain the identity 



(l-t p ) m * = det(l-D(hp)tp) , (33) 

valid cycle by cycle in the Euler products (§) for det(l — C). Here "det" refers to 
the [g x g] matrix representation D(hp); as we shall see, this determinant can be 
evaluated in terms of standard characters, and no explicit representation of D{hp) 
is needed. Finally, if a cycle p is invariant under the symmetry subg roup 1~~Lp C G 
of order h p , its weight can be written as a repetition of a fundamental domain cycle 

tp = 4" ( 34 ) 

computed on the irreducible segment that coresponds to a fundamental domain 
cycle. For example, in fig. |3] we see by inspection that t±2 = and = tf . 

We conclude this section with a few comments about the role of symmetries in 
actual extraction of cycles. In the example at hand, the iV-disk billiard systems, a 
fundamental domain is a sliver of the iV-disk configuration space delineated by a 
pair of adjoining symmetry axes, with the directions of the momenta indicated by 
arrows. The flow may further be reduced to a return map on a Poincare surface 
of section, on which an appropriate transfer operator may be constructed. While 
in principle any Poincare surface of section will do, a natural choice in the present 
context are crossings of symmetry axes. 

In actual numerical integrations only the last crossing of a symmetry line needs 
to be determined (using for example the method of ref. pH|). The cycle is run in 
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global coordinates and the group elements associated with the crossings of sym- 
metry lines are recorded; integration is terminated when the orbit closes in the 
fundamental domain. Periodic orbits with non-trivial symmetry subgroups are par- 



ticularly easy to find since their points lie on crossings of symmetry lines |p9|, j34}] . A 



multi-point-shooting method combined with Newton-Raphson iteration has proven 



very efficient [50] in practice. 



5.1 Boundary orbits 

Before we can turn to a presentation of the factorizations of zeta functions for the 
different symmetries we have to discuss a peculiar effect that arises for orbits that 
run on a symmetry line that borders a fundamental domain |60, 26, 27]. In our 



3-disk example, no such orbits are possible, but they exist in other systems, such 
as in the bounded region of the Henon-Heiles potential and in 1-d maps. For the 
symmetrical 4-disk billiard, there are in principle two kinds of such orbits, one kind 
bouncing back and forth between two diagonally opposed disks and the other kind 
moving along the other axis of reflection symmetry; the latter exists for bounded 
systems only. While there are typically very few boundary orbits, they tend to be 
among the shortest orbits, and their neglect can seriously degrade the convergence 
of cycle expansions, as those are dominated by the shortest cycles. 

While such orbits are invariant under some symmetry operations, their neigh- 
borhoods are not. This affects the stability matrix J p of the linearization perpen- 
dicular to the orbit and thus the eigenvalues. Typically, e.g. if the symmetry is a 
reflection, some eigenvalues of J p change sign. This means that instead of a weight 
1/ det(l — J p ) as for a regular orbit, boundary cycles also pick up contributions 
of form 1/ det(l — hJ p ), where h is a symmetry operation that leaves the orbit 
pointwise invariant; see for example sect. |6.3| . 

Consequences for the zeta function factorizations are that sometimes a boundary 
orbit does not contribute. A derivation of a dynamical zeta function (||) from 
a determinant like ([5|) usually starts with an expansion of the determinants of 
the Jacobian. The leading order terms just contain the product of the expanding 
eigenvalues and lead to the zeta function (|6[). Next to leading order terms contain 
products of expanding and contracting eigenvalues and are sensitive to their signs. 
Clearly, the weights t p in the zeta functions will then be affected by reflections in 
the Poincare surface of section perpendicular to the orbit. In all our applications it 
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was possible to implement these effects by the following simple prescription. 

If an orbit is invariant under a little group TL P = {e,b 2 , ■ ■ ■ ,bh}, then the corre- 
sponding group element in (^) will be replaced by a projector. If the weights are 
insensitive to the signs of the eigenvalues, then this projector is 

1 h 

9p = -rY, h i- ( 35 ) 
n i=i 

In the cases that we have considered, the change of sign may be taken into account 
by defining a sign function e p (g) = ±1, with the "-" sign if the symmetry element 
g flips the neigborhood. Then fl35|) is replaced by 

1 h 

&> = t £ e ( & *) b * ■ ( 36 ) 



The resulting zeta functions agree with the ones given by Lauritzen||26||. We illus- 
trate the above in sect. |6.3| by working out the full factorization for the 1-dimensional 
reflection symmetric maps. 



6 Factorizations of zeta functions 

In the above we have shown that a discrete symmetry induces degeneracies among 
periodic orbits and decomposes periodic orbits into repetitions of irreducible seg- 
ments; this reduction to a fundamental domain furthermore leads to a convenient 
symbolic dynamics compatible with the symmetry, and, most importantly, to a fac- 
torization of zeta functions. This we now develop, first in a general setting and then 
for specific examples. 



6.1 Factorizations of dynamical zeta functions 

According to (|33| ) and (|34]), the contribution of a degenerate class of global cycles 
(cycle p with multiplicity m p = g/h p ) to a zeta function is given by the correspond- 
ing fundamental domain cycle p: 

(l-t^) 9 ^ = det(l-D(h p )t p ) (37) 

Let D{h) = ® a d a D a {h) be the decomposition of the matrix representation D{h) 
into the d a dimensional irreducible representations a of a finite group G. Such 
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decompositions are block-diagonal, so the corresponding contribution to the Euler 
product (H) factorizes as 

det(l - D(h)t) = [] det(l - D a (h)t) da , (38) 

a 

where now the product extends over all distinct d a - dimensional irreducible repre- 
sentations, each contributing d a times. For the cycle expansion purposes, it has 
been convenient to emphasize that the group-theoretic factorization can be effected 
cycle by cycle, as in (|3~?D ; but from the transfer operator point of view, the key 
observation is that the symmetry reduces the transfer operator to a block diagonal 
form; this block diagonalization implies that the zeta functions @ factorize as 

7 = II i > T = Il det (! - D a {hp)tp) . (39) 

S a Set Sa ft 

Determinants of <i-dimensional irreducible representations can be evaluated us- 
ing the expansion of determinants in terms of traces, 

det(l + M) = l + trM+ - ((trM) 2 -trM 2 ) 

+ i ((trM) 3 -3(trM)(tr M 2 ) + 2 trM 3 ) 

+ --- + l((trM) rf ----) , (40) 

and each factor in ( j38|) can be evaluated by looking up the characters Xa(h) = 
tr D a (h) in standard tables[[56|] . In terms of characters, we have for the 1-dimensional 
representations 

det(l - D a {h)t) = 1- Xa (h)t , 

for the 2-dimensional representations 

det(l - D a (h)t) = 1 - Xa {h)t + - ( X a{h) 2 - Xa {h 2 )) t 2 , 

and so forth. These expressions can sometimes be simplified further using standard 
group-theoretical methods. For example, the | ((trM) 2 — trM 2 ) term in QUID is 
the trace of the antisymmetric part of the M x M Kronecker product; if a is a 
2-dimensional representation, this is the A2 antisymmetric representation, so 

2-dim: det(l - D a (h)t) = 1 - Xa(h)t + XA 2 (h)t 2 . (41) 

In the fully symmetric subspace trD^^h) = 1 for all orbits; hence a straight- 
forward fundamental domain computation (with no group theory weights) always 
yields a part of the full spectrum. In practice this is the most interesting subspec- 
trum, as it contains the leading eigenvalue of the transfer operator. 
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6.2 Factorizations of functional determinants 



Factorization of the full functional determinant (f|) proceeds in essentially the same 
manner as the factorization of dynamical zeta functions outlined above. By (KJ3) 
and ( |52]) the trace of the transfer operator C splits into the sum of inequivalent 
irreducible subspace contributions ^ a tr C a , with 



tr C a — d a Xa(h) / dxC(h 1 



X, X) 



This leads by standard manipulations [|i~8| , [46f to the factorization of (|5]) into 



Z z) 



ZJz) 



oo 

eXP, "^Sr"|det (l-Jj) 



(42) 



where Jp = hpjp is the fundamental domain Jacobian. Boundary orbits require 
special treatment, discussed in sect. |5.1| , with examples given in the next section as 
well as in sect. [7|. 



The factorizations (|39|) , (j42|) are the main result of this paper. We now proceed 
to examplify it by a few cases of physical interest. Additional examples (factoriza- 
tion and cycle expansions for the cyclic symmetry groups C3 and Cn) are given in 
ref. f22|. 



6.3 Reflection symmetric 1-d maps: C2 factorization 

Consider /, a map on the interval with reflection symmetry /(— x) = —f(x). Denote 
the reflection operation by Cx = —x. The symmetry of the map implies that if {x n } 
is a trajectory, than also {Cx n } is a trajectory because Cx n+ i = Cf(x n ) = f(Cx n ) . 
The dynamics can be restricted to a fundamental domain, in this case to one half of 
the original interval; every time a trajectory leaves this interval, it can be mapped 
back using C. 

To compute the traces of the symmetrization and antisymmetrization projection 
operators (^8|), we have to distinguish three kinds of cycles[|6l[]: asymmetric cycles 
a, symmetric cycles s built by repeats of irreducible segments s, and boundary cy- 
cles b. The Fredholm determinant can be formally written as the product over the 
three kinds of cycles: det(l — C) = det(l — £) a det(l — £) s det(l — £)&. 
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Asymmetric cycles: A periodic orbits is not symmetric if {x a }n{Cx a } = 0, where 
{x a } is the set of periodic points belonging to the cycle a. Thus C generates a second 
orbit with the same number of points and the same stability properties. Both orbits 
give the same contribution to the first term and no contribution to the second term 
in (28|); as they are degenerate, the prefactor 1/2 cancels. Resumming as in the 



derivation of (|7|), we find that asymmetric orbits yield the same contribution to the 
symmetric and the antisymmetric subspaces: 



det(i-4t) a = TT TT 1-tt > ta 




Symmetric cycles: A cycle s is reflection symmetric if operating with C on the 
set of cycle points reproduces the set. The period of a symmetric cycle is always 
even (n s = 2rig) and the mirror image of the x s cycle point is reached by traversing 
the irreducible segment s of length rig, f ns (x s ) = Cx s . 5(x — f n (xj) picks up 
2rig contributions for every even traversal, n = rrig, r even, and 5(x + f n {x)) for 
every odd traversal, n = rrig, r odd. Absorb the group-theoretic prefactor in the 
stability eigenvalue by defining Ag = —Df ns (x s ), where Df ns (x s ) is the stability 
computed for a segment of length rig. Restricting the integration to the infinitesimal 
neighborhood of the s cycle, we obtain the contribution to trL±: 

z n tr£l -> f dxz n \ (6{x-f n (x)) ±S(x + f n (x))) 

JVs ^ 



n 



/even ±r odd j.r 

/ r b g i f L S 



71 



^ A (±tgy 
~ s h n,ms i-i/A|" 



Substituting all symmetric cycles s into det(l — C±) and resumming as in (0), we 
obtain: 



det(i-/: ± ) s =n 5(^1 



Boundary cycles: In the example at hand there is only one cycle which is neither 
symmetric nor antisymmetric, but lies on the boundary of the fundamental domain, 
the fixed point at the origin. Such a cycle contributes simultaneously to both 
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S(x - f n (x)) and S(x + f n (x)): 

z n tiC n ± -> f dxz n \ (5(x- f n {x))±5(x + f n (x))) 

JVt 2i 

°° r 1 / 1 1 \ 



2 V1-1/AJ 1 + 1/AL 

z n tr£™ -> 



00 00 1 j.r 

6 . _n j.„ . r 6 _ 

2r • 
b 



feS- 1/Ar ' - r t! ™A£ 1 - 1/A| 



Boundary orbit contributions to the factorized Fredholm determinants follow by 
resummation: 

det(l - C + ) b = fi f 1 " 4) , det(l - = fi f 1 - J^r) 

fc=0 V A 6 / fe=0 V A 6 / 

Only even derivatives contribute to the symmetric subspace (and odd to the anti- 
symmetric subspace) because the orbit lies on the boundary. The symmetry reduced 
Fredholm determinants follows by collecting the above results: 

M>)=nn(i-p]nfi(i-p)ii(i-p) 

a k=0\ A a/ S k=0\ A s/fc=0\ A 6 / 

^)=nfi(i-s)nfi(i+|)n(i-Tsr) m 

a k=Q \ iS W S k=0 \ 1V s / k=0 \ 1V b / 

As reflection symmetry is essentially the only discrete symmetry that a map of 
the interval can have, this example completes the group-theoretic factorization of 
determinants and zeta functions for 1-dimensional maps. A specific example is 
worked out in Ref. [BUI 



7 Examples of symmetry induced factorizations 

We conclude with several explicit examples of group theory factorizations of cycle 
expansions of dynamical zeta functions. These expansions are a prerequisite for ap- 
plications of periodic orbit theory to the evaluation of classical and quantal spectra; 
in particular, they were used in the calculations of refs. [§, Q . 
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7.1 Ci factorizations 

As the simplest example of implementing the above scheme consider the C% sym- 
metry which arises, for example, in the Lorenz system p5], in the 3-dimensional 



anisotropic Kepler problem |24j, |30| , pl| or in the cycle expansions treatments of the 
Ising model^]. For our purposes, all that we need to know here is that each orbit 
or configuration is uniquely labelled by an infinite string {s^}, Sj = +, — and that 
the dynamics is invariant under the + <-> — interchange, i.e., it is C2 symmetric. 
In the Lorenz system, the labels + and — stand for the left or the right lobe of the 
attractor and the symmetry is a rotation by n around the z-axis. Similarly, the Ising 
Hamiltonian (in the absence of an external magnetic field) is invariant under spin 
flip. The C2 symmetry cycles separate into two classes, the self-dual configurations 

H — , + H , + + H , H 1 h H — , ■ ■ •, with multiplicity m p = 1, and 

the asymmetric configurations +, — , + H — , h, ■ - •, with multiplicity m p = 2. 

For example, as there is no absolute distinction between the "up" and the "down" 

spins, or the "left" or the "right" lobe, t + = t-, £++_ = , and so on. 

The symmetry reduced labelling p; e {0, 1} is related to the standard Si G 
{+, — } Ising spin labelling by 

If Si = s,_i then Pi = 1 

If Si ^ s^ then p { = (44) 
For example, + = ••• + + + H — • maps into • • • 111 • • • = 1 (and so does — ), 



— h = 1 1 maps into • • - 000 • • • = 0, - + H — = hH hH 

maps into • • • 0101 • • • = 01, and so forth. A list of such reductions is given in table 
4. 

Depending on the maximal symmetry group 7i p that leaves an orbit p invariant 
(cf. sects. H, U), the contributions to the zeta function factor as 

H p = {e}: (l-t p ) 2 = (1-^(1-*-) 
H p ={e,a}: (1-tj) = (1 - t p )(l + t p ) , (45) 

For example: 

^++- = {e}: (l-t++-) 2 = (l-*ooi)(l-*ooi) 

H+- = {e,a}: (1 - = (1 - t ) (l + t ), t+- = t\ 
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This yields two binary expansions. The A± subspace zeta function is given by the 
standard binary expansion (||). The antisymmetric A 2 subspace zeta function (a 2 
differs from only by a minus sign for cycles with an odd number of O's: 

l/C M = (l + to)(l-*i)(l+*io)(l-*ioo)(l+*ioi)(l + «iooo) 
(1 — £1001) (1 + £1011) (1 — £10000) (1 + £10001) 
(1 + tiooio)(l — £ioon)(l — £ioioi)(l + £10111) • • • 
= 1 + to — £1 + (£10 — £i£o) — (£100 — £io£o) + (£101 — £io£i) 

— (£1001 ~~ £i£ooi ~~ £ioi£o + £io£o£i) ~~ (46) 

Note that the group theory factors do not destroy the curvature corrections (the 
cycles and pseudo cycles are still arranged into shadowing combinations). 

If the system under consideration has a boundary orbit (cf. sect. [5.1|) with 
group-theoretic factor h p = (e + er)/2, the boundary orbit does not contribute to 
the antisymmetric subspace 

A 1 A 2 

boundary: (1 - t p ) = (1 - t p )(l - 0t p ) (47) 



This is the l/( part of the boundary orbit factorization of sect. B73 



7.2 3-disc pinball: C$ v factorization 

The next example, the C% v symmetry, can be worked out by a glance at fig. [3]a. 
For the symmetric 3-disk pinball the fundamental domain is bounded by a disk 
segment and the two adjacent sections of the symmetry axes that act as mirrors (see 
fig. |]b). The three symmetry axes divide the space into six copies of the fundamental 
domain. Any trajectory on the full space can be pieced together from bounces in 
the fundamental domain, with symmetry axes replaced by flat mirror reflections. 
The binary {0, 1} reduction of the ternary three disk {1,2,3} labels has a simple 
geometric interpretation: a collision of type reflects the projectile to the disk it 
comes from (back-scatter), whereas after a collision of type 1 projectile continues 
to the third disk. For example, 23 = • • • 232323 • • • maps into • • • 000 • • • = (and 
so do 12 and 13), 123 = • • • 12312 • • • maps into • • • 111 • • • = 1 (and so does 132), 
and so forth. A list of such reductions for short cycles is given in table 5. 
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C^v has two one-dimensional irreducible representations, symmetric and anti- 
symmetric under reflections, denoted A\ and A 2 , and a pair of degenerate two- 
dimensional representations of mixed symmetry, denoted E. The contribution of 
an orbit with symmetry g to the l/£ Euler product (|38| ) factorizes according to 

det(l -D(h)t) = (l- XAi (h)t)(l- X A 2 (h)t)(l-XE(h)t + XA 2 (h)t 2 ) 2 m 

with the three factors contributing to the C% v irreducible representations A\, A 2 
and E, respectively, and the 3-disk zeta function factorizes into ( = CAiObCl- 



Substituting the C 3v characters J56| 





A 1 A 2 E 


e 

r< r~i2 

a v 


1 1 2 
1 1 -1 
1 -1 



into (£H), we obtain for the three classes of possible orbit symmetries (indicated in 
the first column) 



e 

3, ^3 



c,,c. 



A, A 2 E 
'l-t p )(l-t p )(l-2t p + t 
[l-t p )(l-t p )(l + t p + tj 
[l- tp )(l + t p )(l + 0t p -tl) 



2\2 
p) 

2\2 
p) 

2\2 
p) 



(49) 



where a v stands for any one of the three reflections. 

The Euler product @ on each irreducible subspace follows from the factoriza- 
tion fl4~9|). On the symmetric A\ subspace the (a 1 is given by the standard binary 
curvature expansion (||). The antisymmetric A 2 subspace (a 2 differs from only 
by a minus sign for cycles with an odd number of O's, and is given in (f46"l). For the 
mixed-symmetry subspace E the curvature expansion is given by 



I/O 



1 + zt x + z%){1 - z%){\ + zH W0 + z b ti 00 )(l - zH 



10 J 



1 + z 4 tiooi + ^ 8 ^iooi)(l + ^ioooo + z 10 t 2 



'iooi; 

(1 + z 5 tioioi + z txoiox)0- — ^ 5 ^ioon) 2 
1 + zt x + z 2 (t? - tg) + z 3 (t oi - titjj) 



ioooo ; 



+z q 



+z a 



'01 



^0011 + (^001 ~~ tit Q )ti — t\ 

^ooooi + ^oioii — 2ioom + (^ooii — ^01)^1 + ~~ *o)^ioo 



+ • {50) 
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We have reinserted the powers of z in order to group together cycles and pseudo- 
cycles of the same length. Note that the factorized cycle expansions retain the 
curvature form; long cycles are still shadowed by (somewhat less obvious) combi- 
nations of pseudocycles. 

Refering back to the topological polynomial (^3j ) obtained by setting t p — 1, we 
see that its factorization is a consequence of the C^ v factorization of the ( function: 

l/( Al = 1 - 2z , 1/Ca 2 = 1 , 1/Ce = l + z, 

as obtained from flD, (|46|) and ([50| ) for t p — 1. 

An example of a system with C$ v symmetry is provided by the motion of a 
particle in the Henon-Heiles potential p2| 



V(r,0) = -r 2 + -r 3 sin(3fl) . 
2 3 

Our coding is not directly applicable to this system because of the existence of 
elliptic islands and because the three orbits that run along the symmetry axis cannot 
be labelled in our code. However, since these orbits run along the boundary of the 
fundamental domain, they require the special treatment |26j discussed in sect. f>A\ 



Their symmetry is K = {e, a}, so according to (|35f) , they pick up the group- 
theoretic factor hp = (e + er)/2. If there is no sign change in t p , then evaluation of 
det(l - s±zt$) yields 

A l A 2 E 

boundary: (1 - t p f = (1 - t p )(l - Ot p )(l - t p f , t p = t p . (51) 

However, if the cycle weight changes sign under reflection, t ap = —tp, the boundary 
orbit does not contribute to the subspace symmetric under reflection across the 
orbit; 

A x A 2 E 

boundary: (1 - t p ) 3 = (1 - Ot p )(l - t p )(l - t p ) 2 , t p = t p . (52) 



7.3 C\ v factorization 

If an iV-disk arrangement has symmetry, and the disk visitation sequence is 
given by disk labels {eie2£3 . . .}, only the relative increments pi = ej + i — £j mod iV 
matter. Symmetries under reflections across axes increase the group to C^v and 
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add relations between symbols: {e^} and {N — e«} differ only by a reflection |40|, |48[ . 
As a consequence of this reflection increments become decrements until the next re- 
flection and vice versa. Consider four equal disks placed on the vertices of a square 
(fig||). The symmetry group consists of the identity e, the two reflections cr x , a y 
across x, y axes, the two diagonal reflections 013, 024, and the three rotations C4, 
C 2 and Cf by angles tt/2, 71 and 37r/2. We start by exploiting the C 4 subgroup sym- 
metry in order to replace the absolute labels 6j G {1,2, 3, 4} by relative increments 
Pi G {1,2,3}. By reflection across diagonals, an increment by 3 is equivalent to an 
increment by 1 and a reflection; this new symbol will be called 1. Our convention 
will be to first perform the increment and then to change the orientation due to the 
reflection. As an example, consider the fundamental domain cycle 112. Taking the 
disk 1 — > disk 2 segment as the starting segment, this symbol string is mapped into 



the disk visitation sequence l + i2 +1 3 +2 l . . . = 123, where the subscript indicates the 
increments (or decrements) between neighboring symbols; the period of the cycle 



112 is thus 3 in both the fundamental domain and the full space. Similarly, the 



cycle 112 will be mapped into l + i2_il_ 2 3_i2 + i3 + 2l = 121323 (note that the fun- 
damental domain symbol 1 corresponds to a flip in orientation after the second and 
fifth symbols); this time the period in the full space is twice that of the fundamen- 
tal domain. In particular, the fundamental domain fixed points correspond to the 
following 4-disk cycles: 

4-disk reduced 

12 <-> 1 
1234 «-> 1 

13 <-> 2 

Conversions for all periodic orbits of reduced symbol period less than 5 are listed 
in table 6. 



While there is a variety of labelling conventions |57|, [48] , |39fl for the reduced 
C^v dynamics, we prefer the one introduced here because of its close relation to 
the group-theoretic structure of the dynamics: the global 4-disk trajectory can be 
generated by mapping the fundamental domain trajectories onto the full 4-disk 
space by the accumulated product of the C iv group elements g% = C, g 2 = C 2 , 
gi = UdiagC = cr ax i S , where C is a rotation by tt/2. In the 112 example worked 
out above, this yields giu = gig\g\ = C 2 Ca ax i S = adiag, listed in the last column 
of table 6. Our convetion is to multiply group elements in the reverse order with 



28 



respct to the symbol sequence. We need these group elements for our next step, 
the zeta function factorizations. 

The C± v group has four one- dimensional representations, either symmetric (A\) 
or antisymmetric (A 2 ) under both types of reflections, or symmetric under one and 
antisymmetric under the other (B\, B 2 ), and a degenerate pair of two-dimensional 
representations E. Substituting the C iv characters 





A l 


A 2 


B x 


B 2 


E 


e 


1 


1 


1 


1 


2 


c 2 


1 


1 


1 


1 


-2 


C*4, C| 


1 


1 


-1 


-1 





®axes 


1 


-1 


1 


-1 





&diag 


1 


-1 


-1 


1 






into (^) we obtain: 



hp 








A, 


A 2 


Br 


B 2 




E 


e 


(1 


tp, 


8 = 


(i 


-tp) 


(1 - tp) 


(i-*fO 


a-tp) 


(1 


- tpY 


c 2 


(1 




4 = 


(i 


-tp) 


(1 " tp) 


(1-tp) 


(1 " tp) 


(1 


+ tpY 


C*4, C| 


(1 


-€ 


2 = 


(i 


-tp) 


(1 " tp) 


(1+g 


(l + t p ) 


(1 


+ tj) 2 


@axes 


(1 


-« 


4 = 


(i 


-tp) 


(l + t p ) 


(1-tp) 


(1 + tp) 


(1 


-4) 2 


@diag 


(1 


-« 


4 = 


(i 


-tp) 


(1 + tp) 


(1+tp) 


(1-tp) 


(1 


- tl? 



The possible irreducible segment group elements hp are listed in the first column; 
&axes denotes a reflection across either the x-axis or the y-axis, and <Jd%ag denotes a 
reflection across a diagonal (see fig. In addition, degenerate pairs of boundary 
orbits can run along the symmetry lines in the full space, with the fundamental 
domain group theory weights h p = (C 2 + a x )/2 (axes) and h p = (C 2 + cri 3 )/2 
(diagonals) respectively: 

A Y A 2 B x B 2 E 
axes: (1 - t 2 -) 2 = (1 - t p )(l - Ot p )(l - t p )(l - Ot p )(l + t p ) 2 
diagonals: (1 — t|) 2 = (1 - t p )(l - Ot p )(l - Ot p )(l - t p )(l + t p ) 2 (53) 

(we have assumed that t p does not change sign under reflections across symmetry 
axes). For the 4-disk arrangement considered here only the diagonal orbits 13, 24 
occur; they correspond to the 2 fixed point in the fundamental domain. 
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The A\ subspace in C± v cycle expansion is given by 

l/U = (l-to)(l-*l)(l-*2)(l-*Ol)(l-*02)(l-<12) 

(1 - t 001 )(l - t 002 )(l - *01l)(l - *012)(1 " *02l)(l " *022)(1 " *11 2 ) 
(1 — tl22j(l — ^000l)(l — ^0002j(l — £o01l)(l — ^0012)(1 — ^002l) • • • 
— 1 — to — ti — ti — (£q1 — ^O^l) — (^02 — ^0^) — (^12 — M2) 

— (£ooi — ^01) — (^002 — ^02) — (ton — ti^oi) 

— (^022 — ^2^02) — (^112 — t\t\2) — (^122 — ^2^12) 

— (^012 + ^021 + ^0^2 — ^0^12 — ^1^02 — ^Ol) • • • (54) 

(for typographical convenience, 1 is replaced by in the remainder of this section). 
For one- dimensional representations, the characters can be read off |63| the symbol 
strings: Xa 2 {K) = (-1)"°, XbAK) = Xb 2 {K) = (-l) no+n \ where n and 

n\ are the number of times symbols 0, 1 appear in the p symbol string. For B2 all 
t p with an odd total number of 0's and l's change sign: 

1/Cb 2 = (l+to)(l + ti)(l-*2)(l-*oi)(l+fo2)(l + *ia) 

(1 + t 001 )(l - t 002 )(l + toii)(l - *oi2)(l - *<m)(l + *02 2 )(1 " tm) 
(1 + ti2 2 )(l ~ £0001) (1 + ^0002) (1 — ioou)(l + ^0012) (1 + ^002l) • • • 

= 1 + t + t X - t 2 - (t l - t ti) + (t 02 - t t 2 ) + (ti2 - tit 2 ) 
+ (^001 — ^0^0l) — (^002 — ^0^02) + (^011 — tltoi) 
+ (^022 — ^2^02) — (^112 — ^1^12) + (^122 — ^2^12) 

— {toi2 + ^021 + tot\t2 — tot\2 — tito2 ~ ^Ol) • • • (55) 

The form of the remaining cycle expansions depends crucially on the special role 
played by the boundary orbits: by (Q) the orbit £ 2 does not contribute to A 2 and 
Bi, 

1/Ca 2 = (l + to)(l-ti)(l+toi)(l + to 2 )(l-ti2) 

(1 - tooi)(l - i 002 )(l + *01l)(l + t 0V2 )(l + t 2l)(l + *022)(1 " Hl2) 
(1 — ti2 2 )(l + ^000l)(l + ^0002)(1 — io01l)(l — ^0012j(l ~~ ^002l) • • • 
= 1 + to — tl + (^01 — toti) + io2 — t\2 

— (^ooi — ^o^oi) ~~ (^002 — ^0^02) + (^011 + ^l^Ol) 

+^022 — t\22 — {t\!2 — titl2) + (£oi2 + ^021 — ^0^12 — M02) • • • (56) 
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and 



1/Cbx = (1 - t ) (l + *i)(l+*oi) (l-W (l + tia) 

(1 + t 001 )(l - t 002 )(l - *on)(l + *oia)(l + *02i)(l - *02 2 )(1 - *112) 

(1 + t 122 )(l +toOOl)(l - *0002)(1 - *001l)(l + *001 2 )(1 + *002l) • • • 
— 1 — to + t 1 + (toi — toti) — to2 + ti2 

+ (^001 — totoi) — (^002 — ^0^02) — (^011 — titoi) 

—t022 + £l22 — (^112 — £1^12) + (^012 + ^021 — tot 12 — ^1^02) • • • (57) 

In the above we have assumed that t 2 does not change sign under C^ v reflections. 
For the mixed-symmetry subspace E the curvature expansion is given by 

I/Ce = l + h + (-t 2 + h 2 ) + (2i 002 - W - 2t 112 + t 2 tx 2 ) 



+ (2toon — 2t o22 + 2t 2 too2 — ^oi 2 — t 2 2 + 2t 1122 — 2t 2 t 112 

2 + 2+ 2\ 



+ti2 —tot\ ) + (2t 0002 — 2t 0112 + 2t 2 to011 — 2t 0121 — 2t 0211 
+2^00222 — 2t 2 to022 + 2^01012 + 2^01021 _ 2toll02 — ^Ol 2 + 2^02022 
— t 2 to2 2 + 2tim 2 — 2tn 222 + 2t 2 tn 22 — 2ti2122 + ^2^12 2 — t 2 to 2 tl 2 

+2t o 2 (-to 2 + ti 2 ) - 2t 112 (-t 2 + ti 2 )) (58) 

A quick test of the ( = (a 1 (a 2 (b 1 (b 2 (e factorization is afforded by the topolog- 
ical polynomial; substituting t p = z np into the expansion yields 

l/( Al = 1 - 3z , 1/Ca 2 = 1Kb, = 1 , 1/(b 2 = 1/Ce = 1 + z, 
in agreement with (|23|) . 



7.4 factorization 

An arrangement of four identical disks on the vertices of a rectangle has C 2tI symme- 
try (fig.||). C 2v consists of {e, a x , a y , C 2 }, i.e., the reflections across the symmetry 
axes and a rotation by n. C 2 „ is the symmetry of several systems studied in the lit- 
erature, such as the stadium biHiard[[S|, and the 2-dimensional anisotropic Kepler 
potential pi 



This system affords a rather easy visualisation of the conversion of a 4-disk 
dynamics into a fundamental domain symbolic dynamics. An orbit leaving the 
fundamental domain through one of the axis may be folded back by a reflection 
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on that axis; with these symmetry operations go — o~ x and g\ = a y we associate 
labels 1 and 0, respectively. Orbits going to the diagonally opposed disk cross the 
boundaries of the fundamental domain twice; the product of these two reflections 
is just C 2 = o~ x a y , to which we assign the label 2. For example, a ternary string 
0010201. ..is converted into 12143123. . ., and the associated group-theory weight 
is given by . . . gig g 2 g gig go- 

Short ternary cycles and the corresponding 4-disk cycles are listed in table 7. 
Note that already at length three there is a pair of cycles (012 = 143 and 021 = 142) 
related by time reversal, but not by any C 2v symmetries. 

The above is the complete description of the symbolic dynamics for 4 sufficiently 
separated equal disks placed at corners of a rectangle. However, if the fundamental 
domain requires further partitioning, the ternary description is insufficient. For 
example, in the stadium billiard fundamental domain one has to distinguish between 
bounces off the straight and the curved sections of the billiard wall; in that case there 



exists evidence |3| that five symbols suffice for constructing the covering symbolic 
dynamics. 

The group C 2v has four one-dimensional representations, distinguished by their 
behaviour under axis reflections. The A\ representation is symmetric with respect 
to both reflections; the A 2 representation is antisymmetric with respect to both. 
The B\ and B 2 representations are symmetric under one and antisymmetric under 
the other reflection. The character table is 



c 2v 


At 


A 2 


Bi 


B 2 


e 


1 


1 


1 


1 


c 2 


1 


1 


-1 


-1 


O'x 


1 


-1 


1 


-1 




1 


-1 


-1 


1 



Substituted into the factorized determinant fl38f), the contributions of periodic 
orbits split as follows 



9p 








At 


A 2 


Si 


B 2 


e: 


(1 




4 = 


(1 " tp) 


(1 " tp) 


(i-tp) 


(1-tp 


C 2 : 


(1 


"« 


2 = 


(1-tp) 


(l-*p) 


(i-tp) 


{1-tp 


O'x- 


(1 


"« 


2 = 


(1 - tp) 


(1 + tp) 


(i-tp) 


(1 + tp 


O-y. 


(1 


"« 


2 = 


(1 - tp) 


{1 + tp) 


{l + tp) 


(1-tp 
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Cycle expansions follow by substituting cycles and their group theory factors from 
table 7. For A\ all characters are +1, and the corresponding cycle expansion is given 
in flSip. Similarly, the totally antisymmetric subspace factorization A 2 is given by 
(|55l), the B 2 factorization of C± v . For Bi all t p with an odd total number of O's and 
2's change sign: 

i/C Bl = (i + to)(i-h)(i+t 2 )(i+t i)(i-t 02 )(i + t 12 ) 

(1 - tooi)(l + t W2 )(l + t 011 )(l - t 012 )(l - t ai)(l + t 022 )(l + t n2 ) 

(1 — £122) (1 + ^OOOl) (1 — ^0002) (1 — ^OOll) (1 + ^0012)(1 + ^002l) • • • 

= 1 + 1 - h + 1 2 + (t 01 - toh) - (t 02 - t t 2 ) + (t 12 - txt 2 ) 

— (tool — totoi) + (^002 — ^0^02) + (toil — titoi) 
+ (t()22 — ^2^02) + (^112 — ^1^12) — (^122 — t 2 t\ 2 ) 

— (toi2 + fei + totit 2 — totu — tito2 — £2^01) • • • (59) 
For B 2 all t p with an odd total number of l's and 2's change sign: 

I/Cbs = (1 - t ) (l + ti)(l+* 2 )(l+*oi)(l + *02) 

(1 + t 0l)(l + t002)(l - t ll)(l - %2)(1 - *02l)(l - *02 2 )(1 + t n2 ) 

(1 + ^122) (1 + ^oooi)(l + £0002) (1 ~~ £oon)(l — £0012) (1 — £0021) • • • 

= 1 - 1 + *i + h + (toi - toh) + (t 02 - toh) - (t 12 - ht 2 ) 

+ (£001 — to^oi) + (£002 — £0^02) — (ton — titoi) 

— (£022 — £2^02) + (£112 — titi 2 ) + (t 122 — t 2 ti 2 ) 

— (£012 + £021 + totit 2 — t ti 2 — tito2 — £2^01) • • • (60) 

Note that all of the above cycle expansions group long orbits together with their 
pseudo-orbit shadows, so that the shadowing arguments for convergence flTBH still 
apply. 

The topological polynomial factorizes as 

7— = 1 — 32; , — = — = — = 1 + 2, 
consistent with the 4-disk factorization 



8 Summary 

The techniques of this paper have been applied to computations of the 3-disk clas- 
sical and quantum spectra in refs. K| ^2], [50], and to a "Zeeman effect" pinball 
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and the x 2 y 2 potentials in refs. |63|. The main lesson of such calculations is 
that if a dynamical system has a discrete symmetry, the symmetry should be ex- 
ploited; much is gained, both in understanding of the spectra and ease of their 
evaluation. Once this is appreciated, it is hard to conceive of a calculation with- 
out factorization; it would correspond to quantum mechanical calculations without 
wave-function symmetrizations. 

In a larger perspective, the factorizations developed above are special cases of 
a general approach to exploiting the group-theoretic invariances in spectra com- 



putations, such as those used in enumeration of periodic geodesies [^7], f|T| for 
hyperbolic billiards [13] and Selberg zeta functions fT4|. 

Reduction to the fundamental domain simplifies symbolic dynamics and elimi- 
nates symmetry induced degeneracies. While the resummation of the theory from 



the trace sums||12|| to the cycle expansions [p3|l does not reduce the exponential 
growth in number of cycles with the cycle length, in practice only the short orbits 
are used, and for them the labor saving is dramatic. For example, for the 3-disk 
pinball there are 256 periodic points of length 8, but reduction to the fundamental 
domain non-degenerate prime cycles reduces the number of the distinct cycles of 
length 8 to 30. 

In addition, as we demonstrate by explicit calculations in ref. || |42| , cycle 
expansions of the symmetry reduced zeta functions converge dramatically faster 
than the unfactorized zeta functions. One reason is that the unfactorized zeta 
function has many closely spaced zeros and zeros of multiplicity higher than one; 
since the cycle expansion is a polynomial expansion in topological cycle length, 
accomodating such behaviour requires many terms. The zeta functions on separate 
subspaces have more evenly and widely spaced zeros, are smoother, do not have 
symmetry-induced multiple zeros, and fewer cycle expansion terms (short cycle 
truncations) suffice to determine them. Furthermore, the cycles in the fundamental 
domain sample phase space more densely than in the full space. For example, for 
the 3-disk problem, there are 9 distinct (symmetry unrelated) cycles of length 7 or 
less in full space, corresponding to 47 distinct periodic points. In the fundamental 
domain, we have 8 (distinct) periodic orbits up to length 4 and thus 22 different 
periodic points in 1/6-th the phase space, i.e., an increase in density by a factor 3 
with the same numerical effort. 

We emphasize that the symmetry factorization (|49|) of the dynamical zeta func- 
tion is intrinsic to the classical dynamics, and not a special property of quantal 
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spectra (in which context it was used before|23[). Our results are not restricted to 
the Hamiltonian systems, or only to the configuration space symmetries; for exam- 



ple, the discrete symmetry can be a symmetry of the Hamiltonian phase space f25 
In conclusion, the manifold advantages of the symmetry reduced dynamics should 
thus be obvious; full space cycle expansions, such as those included in the appendix, 
are useful only for cross checking purposes. 



9 Appendix 

Here we list the 3- and 4-disk cycle expansions for unfactorized zeta functions. They 
are not recommended for actual computations, as the factorized zeta functions yield 
much better spectra, but they might be useful for cross-checking purposes. 

For the 3-disk pinball (assuming no symmetries between disks) the curvature 
expansion (El) is given by: 



1/C = {l-Z 2 t l2 ){l-Z 2 t l3 ){l-Z 2 t 23 ){l-zH 123 ){l-zH 132 ) 

^12123) " ' ' 



;i - z\ 213 )(l - zH 1232 )(l - z\ 323 )(l - z% 



— 1 — z 2 t\ 2 — z 2 t 23 — z 2 t 3 i — z 3 t± 23 — z 3 ti 32 

— Z 4 [(ti 2 i 3 — ti 2 ti 3 ) + (tl232 — ^12^23) + (^1323 — ^13^23)] 

-z 5 [(t 12123 -t 12 t 123 ) + ---] (61) 

The symmetrically arranged 3-disk pinball cycle expansion of the Euler product 
0|) (see table 2 and fig.|l]) is given by: 

1/C = (1-^ 12 ) 3 (1-^ 123 ) 2 (1-^121 3 ) 3 

(1 — Z 5 ti 2 i 23 ) 6 (l — 2 6 tl21213) 6 (l — £ 6 £l21323) 3 • • • 



1 — 3z t\ 2 — 2z t\ 23 — 3z (tl213 — t 12 ) — Qz (tl2123 — ti 2 t\ 23 ) 

2 ) 
123/ 



— Z G (6ti21213 + 3ti21323 + ^12 — 9^12^1213 — t 2 



-6Z 7 (^1212123 + ^1212313 + ^1213123 + ^12^123 ~~ 3 ^12^12123 — £l23^1213y 
-3Z 8 (2ti2121213 + ^12121313 + 2^12121323 + 2^12123123 



+ 2 ^12123213 + ^12132123 + 3 ^12^1213 + ^12^123 



— 6 ^12^121213 ~~ 3 ^12^121323 ~~ 4 ^123^12123 ~~ ^1213) — " ' ' (62) 

For the symmetriclly arranged 4-disk pinball the symmetry group is C^, of order 
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8. The degenerate cycles can have multiplicities 2, 4 or 8 (see table 3): 



i/C 



(l - zh l2 )\i - z 2 t 13 ) 2 (i - z% 23 ) 8 (i - z\ 213 ) 8 (i - z\ 2U ) 4 

(i - z% 2U ) 2 (i - zH 12A3 )\\ - zH 12123 )\\ - zH 1212A )\i - zH 1213A f 

(1 - ^1214 3 ) 8 (l - ^tl2313) 8 (l " ^1241 3 ) 8 " " " (63) 



and the cycle expansion is given by 



i/C 




where in the coefficient to z 6 the abbreviations S$ and ^4 stand for the sums over 
the weights of the 12 orbits with multiplicity 8 and the 5 orbits of multiplicity 4, 
respectively; the orbits are listed in table 2. 
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Figure captions 

1. The symmetries of three disks on an equilateral triangle. The fundamental 
domain is indicated by the shaded wedge. 

2. Some examples of 3-disk cycles, (a) 12123 and 13132 are mapped into each 
other by <t 2 3, the flip across 1 axis; this orbit has degeneracy 6 under C% v 
symmetries. Similarly (b) 123 and 132 and (c) 1213, 1232 and 1323 are 
degenerate under C 3v . The orbits (d) 121212313 and 121212323 are related 
by time reversal but not by any C 3v symmetry. 

3. The scattering geometry for the disk radius/separation ratio a : R = 1 : 2.5. 
(a) the three disks, with 12, 123 and 121232313 cycles indicated, (b) the 
fundamental domain, ie. a wedge consisting of a section of a disk, two segments 
of symmetry axes acting as straight mirror walls, and an escape gap. The 
above cycles restricted to the fundamental domain are now the two fixed 
points and 1 and the 100 cycle. 

4. The symmetries of four disks on a square. The fundamental domain is indi- 
cated by the shaded wedge. 

5. The symmetries of four disks on a rectangle. The fundamental domain is 
indicated by the shaded wedge. 
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Table 1. List of the 3-disk prime cycles up to length 10. Here n is the cycle length, 
M n the number of prime cycles, N n the number of periodic points and S n the number of 
distinct prime cycles under the C$ v symmetry. Column 3 also indicates the splitting of N n 
into contributions from orbits of lengths that divide n. The prefactors in the fifth column 
indicate the degeneracy m p of the cycle; for example, 3-12 stands for the three prime cycles 
12, 13 and 23 related by 27r/3 rotations. Among symmetry related cycles, a representative 
p which is lexically lowest was chosen. The cycles of length 9 grouped by parenthesis are 
related by time reversal symmetry (but not by any other C% v transformation). 



n 


M n 


N n 


s 


m p -p 


1 













2 


6 


12=6-2 


2 


4-12 + 2-13 


3 


8 


24=8-3 


1 


8-123 


4 


18 


84=6-2+18-4 


4 


8-1213 + 4-1214 + 2-1234 + 4-1243 


5 


48 


240=48-5 


6 


8-(12123 + 12124) + 8-12313 
+ 8-(12134 + 12143) + 8-12413 


6 


116 


732=6-2+8-3+116-6 


17 


8-121213 + 8-121214 + 8-121234 
+ 8-121243 + 8-121313 + 8-121314 
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Table 2. List of the 4-disk prime cycles up to length 8. The meaning of the symbols is 
the same as in Table 1 . Orbits related by time reversal symmetry (but no other symmetry) 
already appear at cycle length 5. List of the cycles of length 7 and 8 has been omitted. 
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Table 3. Number of prime cycles for various alphabets and grammars up to length 
10. The first column gives the cycle length, the second the formula (|ll|) for the number 
of prime cycles for complete A^-symbol dynamics, columns three through five give the 
numbers for N = 2, 3 and 4. 
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Table 4. Correspondence between the C2 symmetry reduced cycles p and the standard 
Ising model periodic configurations p, together with their multiplicities m p . Also listed 
are the two shortest cycles (length 6) related by time reversal, but distinct under Ci- 
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Table 5. C% v correspondence between the binary labelled fundamental domain prime 
cycles p and the full 3-disk ternary {1,2,3} labelled cycles p, together with the C$ v trans- 
formation that maps the end point of the p cycle into the irreducible segment of the p 
cycle. The degeneracy of p cycle is m p = 6n.p/n p . The shortest pair of the fundamental 
domain cycles related by time symmetry are the 6-cycles 001011 and 001101. 
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Table 6. C^ v correspondence between the ternary fundamental domain prime cycles 
p and the full 4-disk {1,2,3,4} labelled cycles p, together with the C^ v transformation 
that maps the end point of the p cycle into an irreducible segment of the p cycle. For 
typographical convenience, the symbol 1 of sect. ^3 has been replaced by 0, so that the 
ternary alphabet is {0,1,2}. The degeneracy of the p cycle is m p = 8np/n p . Orbit 2 is 
the sole boundary orbit, invariant both under a rotation by it and a reflection across a 
diagonal. The two pairs of cycles marked by (a) and (b) are related by time reversal, but 
cannot be mapped into each other by C^ v transformations. 
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Table 7. C 2v correspondence between the ternary {0, 1, 2} fundamental domain prime 
cycles p and the full 4-disk {1,2,3,4} cycles p, together with the C 2v transformation 
that maps the end point of the p cycle into an irreducible segment of the p cycle. The 
degeneracy of the p cycle is m p = 4np/n p . Note that the 012 and 021 cycles are related 
by time reversal, but cannot be mapped into each other by C 2v transformations. The full 
space orbit listed here is generated from the symmetry reduced code by the rules given 
in section 7.4, starting from disk 1. 
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